Role of CuO chains in vortex core structure in YBa2Cu307_5 
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The Bogoliubov-deGennes equations are solved for a proximity model for YBa2Cu307_<5 in a 
magnetic field. The model explicitly includes the effects of the one-dimensional CuO chains, whose 
influence on the vortex core structure is studied. The rapid vortex core contraction as a function of 
field which is seen experimentally at low magnetic fields is naturally explained by the presence of 
the chains. 



I. INTRODUCTION 

As first observed by Golubov and Hartmanni us- 
ing scanning tunneling microscopy (STM), the vortex 
cores in NbSe2 shrink with increasing magnetic field. 
Subsequent muon spin rotation (/iSR) experiments con- 
firmed this to be the case in NbSe^ but also in 
CerW, YNi 2 B 2 G^, LuNi 2 B 2 C 5 , V 3 S£, Vl, Nb 3 Sr4 and 
YBa 2 Cu 3 07_,^^iii^. The STM experiments probe the 
spatial variation of the local density of states, whereas 
/iSR is sensitive to the spatial dependence of the local 
internal magnetic field B(r). The vortex core size is de- 
termined from the //SR measurements by fitting to a the- 
oretical function for B(r) that includes a cutoff function 
F(G, £), where G are the reciprocal lattice vectors and £ 
is the superconducting coherence length. The functional 
form of F(G, £) depends on the spatial dependence of the 
superconducting order parameter A(r) in the core region. 
Since there is no way of knowing exactly what this is in 
a real material, the fitted value of £ reflects differences 
between the theoretical model and the real spatial de- 
pendence of the local field about the vortex cores. Con- 
sequently, £ is generally not the coherence length, but 
rather a measure of the vortex core size. A second defini- 
tion of the vortex core size is the radius ro at which the 
supercurrcnt density \j{r)\ calculated from B(r) reaches 
a maximum. While this definition is robust with respect 
to the assumed model for B(r), there is a contribution 
to the field dependence of ro that comes naturally from 
the overlap of the j(r) profiles of neighboring vortices^. 

Kogan and Zhelezinai^ have proposed a model based 
on weak-coupling BCS theory that explains the field de- 
pendence of the core size in clean high-K superconductors 
as being due to a field-dependent superconducting coher- 
ence length. Their model qualitatively describes the /iSR 
results for CeRu 2 , NbSe 2 , V 3 Si and YNi 2 B 2 C. A field- 
dependent coherence length has also been suggested to 
be the source of the anomalous field-independent flux- line 
lattice form factor observed in small-angle neutron scat- 
tering measurements on CeCoIn^. However, Ichioka 
and Machidai£ have recently argued that this is caused 
by paramagnetic moments due to Zeeman splitting of the 
Fermi surfaces for spin- up and spin-down electrons. 



Within the framework of the microscopic theory, the 
field dependence of the vortex core size can be explained 
without invoking a field-dependent coherence length. So- 
lutions of the quasiclassical Usadel equations for a dirty 
s-wave superconductor-i* 2 -, and solutions of the quasiclas- 
sical Eilenberger equations for clean s-wave and d-wave 
superconductor a 17 ! 18 show that the field dependences of 
the electronic and magnetic structures of the vortex cores 
are coupled. As explained in Refs. [133, the shrink- 
ing of the vortex cores with increasing H occurs due to 
an increased overlap of the wave functions of the quasi- 
particle core states from nearest-neighbor vortices. This 
delocalization of quasiparticlcs, beginning with the more 
spatially extended wave functions of the higher-energy 
core states, increases the slope of A(r) near r = 0, which 
corresponds to a reduction in the size of the vortex core. 
Experimentally, this picture is strongly supported by the 
remarkable correlation found in V 3 Si£ and NbSe^ be- 
tween the field dependences of the core size and the elec- 
tronic thermal conductivity. 

As pointed out in Ref. [l3| the /iSR measurements of 
NbSe 2 and YBa 2 Cu 3 07_,5 are unusual in that at low 
fields the core size £ exceeds the value of the coher- 
ence length calculated from the upper critical field i? c2 . 
Vortex cores larger than estimated from H C 2 have also 
been observed at low field by STM on the 7r-band of 
the two-gap superconductor MgB 2 20 . While supercon- 
ductivity on both the ir and a bands of MgB 2 con- 
tribute to the electronic structure of the vortex cores, 
at low field the dominant contribution comes from the 
loosely bound quasiparticle core states associated with 
the smaller gapped 7r-band 21 . With increasing H these 
core states rapidly delocalize so that at high field the 
core size, and hence iT c2 , is determined by the intrin- 
sic superconductivity on the er-band. Like MgB 2 , there 
is experimental evidence for distinct energy gaps on dif- 
ferent Fermi sheets in NbSe 2 22 i 23 ' 24 . Recently, the ef- 
fects of the Fermi-surface sheet dependent superconduc- 
tivity on the vortex core size became discernible in a low- 
temperature /iSR study of NbSe^. In the same spirit, 
one of us suggested that the large vortex cores at low 
field in YBa 2 Cu 3 07_«5 may be caused by the occurrence 
of superconductivity on the CuO chain bands. 13 
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While all of the cuprate high temperature supercon- 
ductors are based around conducting two-dimensional 
CuC>2 layers, YBa2Cu 3 07_a and YBa 2 Cu40s are unique 
among the cuprates in having an additional type of 
conducting layer, made of one-dimensional CuO chains. 
Band structure calculations^ 5 - and recent photoemission 
experiment o 26 ^ 27 ' 28 ! 29 show that the chains are far from 
half-filling and are therefore unlikely to be strongly cor- 
related, in contrast to the CuC>2 planes. The metallic 
nature of the chains is inferred, primarily, from trans- 
port and a.c. conductivity measurements.— 

The pairing mechanism for chain superconductivity is 
not well established. Penetration depth anisotropy mea- 
surements have also demonstrated that the chains be- 
come superconducting at the same transition tempera- 
ture as the CuC>2 planes^ Given the significant differ- 
ences in band structure between the two, the most nat- 
ural explanation for the single transition temperature is 
that chain superconductivity arises from the proximity 
effect, mediated by single-electron hopping between the 
chains and planes. A generic feature of YBa2Cu30y_,5 
proximity models is the presence of a small pairing en- 
ergy scale associated with chain superconductivity, which 
is in addition to the large energy scale associated with 
pairing in the CuC>2 planes. The small energy scale man- 
ifests itself, for example, as an inflection point in the 
temperature dependence of the superfluid density 32 i 33 
The absence of this feature in microwave experiments on 
YBa2Cu307_a originally appeared to indicate a failure 
of the proximity model,— but has since been shown to 
be consistent with the fact that a fraction S of oxygen 
sites are vacant in the CuO chains f 3 -^ More recently, /zSR 
experiments on YBa2Cu30j_5 have found an inflection 
point ) 12 ' 35 but the clearest evidence for proximity cou- 
pling of the chains comes from recent fj.SH experiments 
on YBa2Cu40g where there is no chain disorder. 36 

The goal of our work is to demonstrate that the 
low-field vortex core contraction in YBa2Cu30y-5 is 
consistent with the YBa2Cu307_<5 proximity model for 
chain superconductivity and is not due to unconven- 
tional mechanisms related to strong correlations (for ex- 
ample, doping-dependent vortex core expansion in under- 
doped La2- x Sr x Cu04 has been attributed to coexisting 
antiferromagnetism— ). This work is part of a broader 
effort to understand how the CuO chains influence var- 
ious electronic properties, motivated first by the possi- 
bility of novel physics associated with having a metallic 
one-dimensional system coupled to a strongly-correlated 
superconductor, and second by a desire to separate the 
effects of chains from physics related to strong correla- 
tions. 

The idea that the different gap energies in multi- 
band superconductors should introduce distinct magnetic 
field scales has been explored in theoretical models for 
MgB 2 ^^J0JL^. and YBa 2 Cu 3 07_ 5 4 3 . In particular, 
it appears that in these materials the magnetic field de- 
pendence of the density of states (DOS) and related prop- 
erties such as the specific heat can be understood if one 
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FIG. 1: Structure of the bilayer model, (a) The model 
consists of a single plane-chain bilayer with on-site en- 
ergies tj,o (j — 1,2) and single-electron hopping matrix 
elements along nearest-neighbor ((,,„) and next-nearest- 
neighbor (tj t „ nn ) bonds as indicated, (b) The model has two 
bands, with the Fermi surfaces as shown. The model parame- 
ters are {ti,o, ii.nn, ii.nnn} = {1, —1,0.45} for the plane layer, 
{tz.O) £2,nn} = {2.4, —2} for the chain layer, and the interlayer 
hopping amplitude is t± — 0.3. The pair interaction in the 
plane layer is V = 1.6, which produces a zero-field pair am- 
plitude A = 0.39 for the single-layer model and A = 0.33 for 
the bilayer model. All energies are in units of |fi, nn |. 



accounts for the presence of both a large and a small su- 
perconducting gap. Furthermore, the idea that the core 
size in different bands should depend on the gap in each 
band has been explored in Refs. HHl], although an ex- 
plicit calculation demonstrating vortex core contraction 
at low fields has not, to our knowledge, been made. 

A brief description of the YBa2Cu307_,5 proxim- 
ity model, followed by derivations of the appropriate 
Bogoliubov-deGennes equations, are presented in Sec. [Hi 
Results of the calculations are given in Sec. lIIIl including 
the main result that the observed vortex core shrinkage 
is indeed consistent with the proximity model for chain 
superconductivity. The results are discussed in a broader 
context in Sec. IIV1 and a brief concluding statement is 
made in Sec. IVl 



II. THEORY 

In this section we derive the Bogoliubov-deGennes 
(BdG) equations appropriate for the proximity model of 
superconductivity in YBa2Cu307_5. Our derivation is 
similar to ones described, for example, in Refs. [44U45| | but 
with the additional complications of multi-band super- 
conductivity. Proximity models for YBa2Cu307_,5 have 
been discussed in detail elsewhere and we refer the reader 
to Refs. [32,34] for more extensive discussions. 

The geometry of the plane-chain model is illustrated 
in Fig. [TJ The model consists of a single bilayer, com- 
prising a two-dimensional layer (aligned with the x-y 
plane) and a layer of one-dimensional chains (aligned 
with the y-axis). The 2D plane represents a Cu02 
layer and is coupled via single-electron hopping to the 
chain layer. This is the simplest model that contains 
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the essential physics of multi-band superconductivity in 
YBa 2 Cu 3 07_,5. There is an intrinsic pairing interaction 
V in the plane, but the chains are intrinsically normal, 
which means that the superconducting order parameter is 
nonzero in the plane layer only. The chains are, nonethe- 
less, superconducting and exhibit a gap in their DOS. 
One important feature of this model that distinguishes 
YBa2Cu307_,5 from other multiband superconductors is 
that the hybridization of the plane and chain layers is 
strongly fc-dependent and, consequently, the induced gap 
in the chain layer does not have a simple <i-wave symme- 
try, even though the order parameter does i 34 ' 46 A con- 
sequence of this is that the chain-projected DOS ex- 
hibits two pairs of coherence peaks , 32 i 34 i 43 i 47 meaning 
that there is more than one superconducting energy scale 
in the chains. We refer to these as the "small" and 
"large" energy scales. 

We consider only magnetic fields aligned with the crys- 
talline c-axis (which we align with the z-axis), perpen- 
dicular to the Cu02 plane, such that screening currents 
circulate within the plane and chain layers. It is for this 
configuration that the vortex core contraction is seen in 
/iSR experiments^ The total Hamiltonian for the model 
can be broken into pieces 



the hopping matrix elements in Eq. ([3]) modified accord- 
ingly. However, our calculations are performed for large 
fields where, as we show below, JB(r) is small and can 
be neglected. 

We take a second-nearest neigbor model with zero-field 
matrix elements tin = t\ t a, tiu,j) = £i,nn, and t-y/uj\\ — 
^i,nnn, where and refer to nearest and next- 

nearest neighbors respectively (cf. Fig. [1]). In the zero- 
field limit, the dispersion for the plane layer is ei(k) = 
t% t o + 2ti jim (cos k x + cos k y ) + 4ii innn cos k x cos k y . 

The local superconducting order parameter A^- is de- 
termined self-consistently under the assumption that the 
pair interaction V is attractive for nearest neighbor elec- 
trons but vanishes otherwise. Then. 



V 



(cufoOerrOi) +c u (r i )ci T (r j ))%,j>. (5) 



The d-wave component, defined by 



A(r,) = A, 



(6) 



is the dominant component of the order parameter. 
The isolated chain layer is described by a Hamiltonian 



H = H 1 +H 2 + H A 



(1) 



H 2 =53* 2 « C L( r i) C 2a(ri) 



(7) 



where Hi is the Hamiltonian for the isolated plane, H^ 
the Hamiltonian for the isolated chains, and H± the 
single-electron hopping term that couples the two layers. 
For comparison, we also consider a single-layer model de- 
scribed by Hi alone. 

The BdG Hamiltonian for the isolated plane is 



Hi = Xl* 1 w c L( r <) Clff ( r i) + X)[ A< J c iT( r <) c u( r J 



Ay c u( r j) c iT( r i)]> 



(2) 



where ci cr (ri) is the annihilation operator for an elec- 
tron in the plane on site i with spin a, and position 
Yi = (xi,yi), tuj are hopping matrix elements, and Ay- 
are superconducting pair energies. The subscripts "1" 
and "2" refer to the plane and chain layers respectively. 
The hopping matrix element tuj between sites i and j 
includes the effects of the magnetic field via the Peierls 
substitution: 



tuj C 



-i(e/Hc) J r r > dr-A(r) 

r. vi+vj , ,1 



(3) 



where tuj are the zero-field matrix elements. Here 
A(rj) = —BoyiX is the static magnetic vector potential, 
where Bq is the uniform applied magnetic field and 



a = eBo/hc. 



(4) 



In principle, the inhomogeneous magnetic field B(r) = 
BqZ + (5B(r) should be calculated self-consistently and 



where tin — ^2,0 and t<nj = t2,nn f° r i and j nearest- 
neighbor sites belonging to the same chain. Note that, 
because of our choice of gauge, the hopping matrix 
elements are unchanged by the magnetic field. The 
zero-field dispersion for the chains is t2(k) = £2,0 + 
nn cos ky. The layers are coupled by interlayer hop- 
ping: 

H± = t± X][ C L( r ;)e2<r( r i) + 4A r i)cia(ri)], (8) 

which mixes the chain and plane wavefunctions. 
Rather than attempt a quantitative description of 
YBa2Cu307_<5, we choose band parameters (cf. Fig. [1]) 
which are optimal for numerical calculations, but which 
preserve the general features of the YBa2Cu307_5 Fermi 
surface. 

While the Hamiltonian is not periodic, there is 
nonetheless a quasi-periodicity which allows us to define 
an L x x L y magnetic supercell containing N = L x L y /a^ 
atomic lattice sites (do is the lattice constant) and en- 
closing an even number of flux quanta, where the super- 
conducting flux quantum is $o = hc/2e. We take two 
vortices per supercell so that 



Bn = 



_2$o_ 

L X Ly 



(9) 



Assuming there are Nk = Nk x Nk y supercells in the sys- 
tem, we can define Bloch states via the transformation 



CniKcr = ^ cw(rj + Rf) 
1=1 



-i(K-Rj+Q!XiYi) 



(10) 
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where Rj = (X], Yj) are the supercell lattice vectors la- 
belled by /, and where r 4 = (xi,yi) now, and hereafter, 
labels sites with site index i £ [1, AT] within the magnetic 
supercell. The supercell wavevector is K = 2Tr(n x , n y )/L, 
where L = Nk x L x = N^ y L y is the linear dimension 
of the system and n x and n y are integers such that 



n x e [1, N k 



£ [l,Nky\- The Hamiltonian is block- 



diagonal in this new basis, and has the form 
H = ^^*J (K)ffj . (K) |. (K) 



(11) 



K ij 



with 



Hij(K) = 



t uj (K) Aij (K) 

A].(K) -t iy (-K) 

t± 

-t± 



t± 



2tf(K) 







-id 





-t 



where wj(K) 



«i(-K)' 
(12) 

and 



^nij (-^-) — ^tj 



(13) 

The hopping matrix elements £ n y(K) have periodic 
boundary conditions at the edges of the supercell: an 
electron at Yj which leaves the supercell via one of its 
edges is periodically mapped back onto site Vi belong- 
ing to the supercell via the appropriate supercell lattice 
vector R = (X,Y). 

In terms of Bloch states, the superconducting gap am- 
plitude is 



A, 



( K ) = ^E^'( K - K ') 



K' 



'(cij-K'iCliK'T +CmK'JC1j-K't)> (14) 



with Vij(q) — — 3^{i,j> e ' q R ' I n many experiments, in 
particular /iSR, it is not the order parameter but the 
magnetic field profile which is measured near a vortex 
core. This is directly related to the current densities in 
the chain and plane layers. The 2D current density at 
a site i in layer n is defined by averaging the current 
densities flowing into and away from the site, 



JS°(ri) = 



2ha, 



-Im < Sr ^( K )( C ^Ka C «jK ( r)(15) 



where the prefactor | comes from the average and Sr = 

Yi - Yj + R. 

We remark that the calculations described above are 
gauge invariant provided Nk x = L y /ao and Nk y = L x /a . 
In practice, it is not feasible to sum over such a large num- 
ber of k-points at low fields where L x and L y are large, 
and by necessity we use a reduced set at the lowest field 
strengths. A consequence of this approximation is that 
the current in the normal state does not vanish identi- 
cally. This is particularly problematic for the chain layer. 
For system sizes up to L x ^ y = 14ap no approximation is 
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FIG. 2: (color online) Self-consistent solutions for the vor- 
tex structure (diamond lattice) with applied field Bo = 
2<E> /2500ag. The d-wave order parameter is shown as a func- 
tion of position (in units of ao) for (a) the single-layer and (b) 
bilayer models. The current magnitudes are shown for (c) the 
single-layer and (d-f) the bilayer models. For the bilayer the 
plane (d), chain (e), and interlayer (f) current amplitudes are 
shown. 



made, while for systems up to L x 



L VtX /2a{j. For the largest system sizes, L x 



40a , Nkx,ky 



50ao and 



L x , y — 60ao, we have taken Nk x ,ky = 5. We have checked 
that the spurious normal-state current in the largest sys- 
tems is at least an order of magnitude smaller than the 
currents reported here in the superconducting state. 



III. RESULTS 

In this section we describe the results of self-consistent 
solutions of the BdG equations for the vortex lattice. Our 
goal is to explain the observed magnetic field dependence 
of the core size at low fields. We have performed calcula- 
tions for diamond, square, and monoclinic lattice struc- 
tures and have found the same qualitative results for the 
vortex core size in all cases. We present results for the 
diamond lattice, for which L x = L y = L and the lattice 
vectors for the primitive unit cell are (L/2, ±L/2). 

The self-consistently determined d-wave gap and cur- 
rent distributions are shown in Fig.[5]near a single vortex. 
Note that the vortex cores shown in Figs. H£a,b) have a 
radius of roughly 2ao, whereas the coherence length is 
~ 5ao in optimally-doped YBCO6.95 ■— This discrepancy 
results from our having taken the order parameter to 
be twice what is appropriate for quantitative models of 
YBa 2 Cu307_«5. We have done this so that energy scale 
of the the induced gap in the chain layer lies near the 
middle of the range of numerically accessible magnetic 
fields. 

Figure[2]illustrates the various effects of proximity cou- 
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pling on the vortex structure. First, there is an over- 
all reduction of the order parameter in the plane layer 
owing to the presence of the chains. This is a general 
feature of proximity models which is independent of the 
magnetic field: while the plane induces superconductiv- 
ity in the chains, the (intrinsically-normal) chains also 
degrade superconductivity in the plane. Similar physics 
has been found in multiband models for MgB2 where im- 
purity scattering between the 7r-band and cr-band mixes 
the two bands42 In our model, this mixing comes from 
the interlayer hopping and means that the vortex core 
in the bilayer model is slightly larger than for the single 
layer model. 

Second, we note that there is an anisotropic suppres- 
sion of the order parameter near the vortex cores which 
is evident in the bilayer model, with the cores being ex- 
tended along the chain direction. This follows from the 
anisotropy of the current in the plane, Fig. [U(d), which 
itself follows from two features of the proximity model: 

(i) the chains and plane carry currents in parallel, and 

(ii) the chains only carry currents in the y direction. The 
current in the plane is consequently larger at positions 
along the y-axis, where it flows entirely in the x direc- 
tion [and the chain-current therefore vanishes, Fig.[2fe)], 
than at corresponding positions along the x-axis. 

Third, there is an interlayer current, Fig. [3Jf ) , that 
has a quadropolar structure and introduces a small 
quadropolar in-plane component to the magnetic field. 
Since the interlayer current density is an order of mag- 
nitude smaller than the intralayer current density, the 
in-plane component is small. 

We wish to extract a vortex core size from our calcula- 
tions; however, the vortex core is not a well-defined ob- 
ject and the core size is not uniquely defined. In simple 
BCS superconductors, the various definitions give sim- 
ilar results^ but the situation is more complicated in 
YBa2Cu307_a where there is more than a single length 
scale. The most obvious measure of the vortex core size 
is the length scale over which the order parameter ap- 
proaches its asymptotic value, usually determined from 
the gradient of the order parameter near the vortex core 
center. 1 - This definition does not work well when the co- 
herence length and the lattice constant are comparable, 
as we have here. Similarly, the commonly-used defini- 
tion that the core size is given by the radius at which the 
current density is a maximum suffers from poor resolu- 
tion due to the discreteness of the atomic lattice. In the 
following, we discuss a procedure for extracting the core 
size from the vorticity of the current distribution. 

This approach is motivated by the fact that /iSR ex- 
periments measure the distribution of the magnetic field 
in the vortex lattice. The magnetic field inhomogene- 
ity <5B(r) = B(r) — Bo can be calculated from the cur- 
rent density profiles via Maxwell's equation V x <5B = 
(47r/c)J, where J(r) is the volume current density. The 
z-component of 5B(r) satisfies 




FIG. 3: (color online) Magnetic field for a single magnetic unit 
cell containing two vortices. Results are shown for the bilayer 
model with So = 24>o/2500aQ . The calculations assume that 
ao = 5 A, d z = 10 A, and that the energy scale is |ii, n n| = 100 
meV. 
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FIG. 4: (color online) Vorticity U)(r) of superfluid currents 
in (a) the plane layer and (b) the chain layer for the bilayer 
model at B = 24>o/576aQ. 



which can be solved using a Jacobian relaxation scheme. 
For illustrative purposes, we make the approximation 
that the small interlayer currents can be neglected and 
that the current is uniformly distributed along the c-axis, 
i.e. that J (r) = [3\ D {x, y) + 3\ D {x 1 y)]/d z is independent 
of z, where d z is the c-axis lattice constant for the atomic 
unit cell. In order to extract quantitative values, we 
adopt approximate model parameters for YBa2Cu307_5 
(see caption of Fig. [3]). The plot of B z (r) shown Fig. [3] is 
for a magnetic field near the lower limit of what is com- 
putationally accessible. We note that the field varies by 
« 0.5% over the magnetic unit cell, which justifies the 
approximation made in Eq. ((3J that the field is uniform. 

One can extract a core size from <5B(r), for exam- 
ple by fitting to a Ginzburg-Landau form for the vor- 
tex lattice^ however, we note from Eq. (fTtj|) that the 
vorticity 



V 2 SB z (r) 



4tt 



w(r) = V x J(r), 



(17) 



(V x J) • z, 



(16) 

gives the vortex core size directly. For the simple ex- 
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ample of an isolated vortex in an isotropic medium, 
J ~ #tanh 2 (r/£)/r, and u>{r) decays exponentially for 
r > £. One can then extract a characteristic vortex core 
size from the second moment of the position along the 
x-axis since 



/ G? 2 rw(r)x 2 
J d 2 r u)(r) 



0.693147C. 



In a vortex lattice, the second moment is not well de- 
fined since w(r) satisfies 



d 2 ruj(r) = 0, 



(18) 



(where the integral is over the vortex unit cell) and is 
therefore not positive definite. However, if the circulation 
around each vortex is positive (negative), then we can 
still extract a characteristic core size based on the region 
over which ui(r) is positive (negative). We define the 
extent of the vortex pa along a direction n as 



PA 



L>o d2ruj ( r )( r - ii 
L>o d2rw ( r ) 



1/2 



(19) 



where the integral is taken over a single vortex. This 
definition is not unique, but it serves to illustrate the 
physics of the vortex core contraction at low fields. The 
main advantage of this approach is that it is relatively 
insensitive to the discreteness of the lattice. 

Examples of the vorticity are shown in Fig. |4] for the 
plane and chain layers. One sees that the vortex core 
in the chain layer is highly anisotropic, and is extended 
along the y-direction. In Fig. EJa), we plot the field- 
dependence of the core sizes pa and p$ for the plane and 
chain layers based on their separate current distributions 
J 2£) (r) and 3\ D (r) respectively. The main point of this 
figure is that while the core size in the plane layer depends 
only weakly on the field, pf (and to a lesser extent pa) 
in the chain layer varies rapidly with Bq for Bq < B*. 
Empirically, 



B* ps 0.001(2$ /ao) 



(20) 



for the model parameters used in this work. In Fig.[S£b), 
we show similar calculations for p x and pf based on the 
combined current J 2D = S\ D + J 2D . Again, there is 
a rapid core contraction with increasing Bq, primarily 
along the y-direction, for Bq < B* . Experimental mea- 
surements of the core sizei^, plotted in Fig. 02c), show 
a similar variation at low field. In comparison, there 
is a relatively weak low-field core-size variation for the 
single-layer model [Fig. [5£b)], a factor of about 1.6 over 
the range of B shown, in quantitative agreement with 
earlier calculations. 17 Figure EJb) is the main result of 
this work. 

We note that the results in Fig. OH are in qualitative 
agreement with a simplified quasiclassical "doppler-shift" 
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FIG. 5: (color online) Core size for the bilayer model along 
the x and y directions for (a) the plane and chain layers and 
(b) for the combined current distribution. The core size for 
the single-layer (SL) model is also shown, (c) Experimental 
data from Ref. [12] is shown for comparison, (d) The chain- 
projected density of states for the bilayer model is plotted for 
different fields (Bq is in units of 2$o/ao)- 



calculation that has been reported previously^. The cur- 
rent results confirm the validity of the previous approxi- 
mate calculations. 

The density of states (DOS) for the chain layer 
[Fig. OHd)] shows that there are two distinct supercon- 
ducting energy scales in the chain spectrum, a large gap 
El ~ 0.35 and a small gap Es ~ 0.1. The two-gap spec- 
trum originates from the one-dimensional structure of the 



chains, and is discussed at length in Refs. [32U34J ] . In 
the range of fields explored (which are much lower than 
the upper critical field), the chain DOS in the interval 
\E\ < Es is a strong function of field for B < B* but sat- 
urates for Bq > B* . This illustrates the close connection 
between Es and B* . We expect that B* is the field at 
which the vortex cores in the chain layer begin to overlap. 
We estimate a BCS length scale £, c hain = fov f, chain /nEs 
for the small gap, where vp, chain is the y-component of 
the Fermi velocity in the chain. For our model parame- 
ters this gives £, c hain ~ lOao, which is close to the low-field 
value of p$ for the chain layer shown in Fig. 03a). For 
the diamond lattice, the vortex spacing along the chain 
direction is \/2£ m , where £ m = \/&o/Bo is the mag- 
netic length. Then, the vortex cores will overlap when 
*/2l m fa i^chain, which gives an estimate for B* of 



B* 



$ 



^cha 



For the current model, this gives B* 
in good agreement with Eq. (j2"D|) . 



(21) 



0.0025(2$ /a§), 
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IV. DISCUSSION 

In this section, we discuss our results in the context of 
related published work. 

A number of tunneling experiments on 
YBa2Cu307_ j 47 i 49 have found a spectrum with multiple 
superconducting energy scales. While the origin of 
the different scales has not been firmly established, 
there is evidence that they arise from a single pairing 
interaction, consistent with the proximity model4£ The 
smallest of the measured gaps is ~ 5 meV in YBCO6.95 
and it is believed to arise from chain superconductiv- 
ity. If we then take Hvf, chain = 4.12 eVA from first 
principles band structure calculations 25 , then we get 
£chain = 262 A. This, using Eq. (f2"Tj) , gives a crossover 
field of B* ~ 1.5 T, which is in remarkably close agree- 
ment with experimental measurements reproduced in 
Pig. \5^c). This has two implications. First, it appears to 
indicate consistency between two distinct experiments, 
one of which (tunneling) is surface sensitive. Second, 
it strengthens the case that the proximity model is 
appropriate for YBa2Cu307_5. 

One of the key features of the YBa2Cu307_a proximity 
model is that the pairing interaction resides within the 
plane layer and that pairing in the chains results from 
single-electron hopping between the physical layers. In 
this model, the induced gap in the chain layer is propor- 
tional to the gap in the planes. A similar model has been 
introduced for MgB2i 39 ' 40 a domininant intraband pair- 
ing interaction in the a-band is assumed, and a subdom- 
inant pairing interaction in the 7r-band arises through 
interband pair-tunneling. There is a qualitative similar- 
ity in the field dependence of the DOS between the two 
models: in both cases, the low-energy DOS fills in rapidly 
as the field increases, but the energy of the gap edge is 
nearly field-independent [Fig. [5(d)]. This should be con- 
trasted with single-band superconductors where the gap 
is field-dependent^ In YBa2Cu307_5, field-induced pah- 
breaking occurs primarily in the chain layer, but the pair- 
ing interaction resides in the plane layer. 

There are also important physical distinctions between 
the MgB2 and YBa2Cu307_5 models. In MgB2 models, 
pairing is generally assumed to occur in the short-ranged 
s-wave channel. Thus, Cooper pairs belong either entirely 
to the cr-band or 7r-band. In contrast, a significant contri- 
bution to chain superconductivity in the YBa2Cu307_,5 
model comes from pairing correlations between electrons 
in the plane and chain layers^ Furthermore, the pre- 
dominantly d-w&ve symmetry of the order parameter in 
the plane layer cannot imply a fourfold-symmetric chain 
gap because the underlying chains are one-dimensional. 
Thus, unlike in MgB2, the k- resolved excitation spectrum 
of the chains is quite complicated and there is no single 
gap energy that one can attach to chain superconductiv- 
ity. 

One consequence of this is that, since B* is associated 
with the smallest of the chain gaps, a significant super- 
fluid density remains in the chain layer when Bq > B* . 



Even at the largest field studied, the chain DOS at the 
Fermi energy is about half its value in the normal state 
[cf. Fig. [SJd)], meaning that Cooper pairs formed with 
binding energies corresponding to the large gap are not 
broken by the magnetic field. 

One interesting question, which is beyond the scope of 
this work, is how the structure of the vortex lattice it- 
self is affected by the presence of the chains. It has been 
found that, at fields less than 4 T, the vortex lattice 
in YBa2Cu307_5 has a distorted hexagonal symmetry 5 * 
with the distortions apparently originating from the CuO 
chains^ At higher fields there is a crossover to a square 
lattice, which is expected from the <i-wave symmetry of 
the order parameter in the Cu02 planes. These experi- 
ments suggest that chain superconductivity is degraded 
for B > 4 T. While this field is approximately 2.5 times 
the value of B* we extracted by eye from the data in 
Fig. [Sfc), we suggest that the crossover in the vortex lat- 
tice structure is generally consistent with both the /iSR 
measurements and the proximity model. Quantitative 
calculations are needed to establish rigorous consistency. 

Finally, we note that while the calculations in this work 
assume that the chains are infinitely long, YBa2Cu307_,5 
has a fraction S of chain-layer oxygen sites which are va- 
cant. O-vacancies effectively divide the chains into frag- 
ments of varying lengths £, and it is worth considering 
how this affects the results presented here. Based on 
our earlier assertion that B* is the field at which vortex 
cores in the chain layer overlap, we suggest that the low- 
field vortex core-size variation should be easily observ- 
able provided that the mean value £ of £ is larger than 
2£,chain ~ 500A. The experiments of Ref. [l^] span a range 
of fillings between YBCO6.57 and YBCO6.95 so it is possi- 
ble that there are large sample-to-sample variations in £. 
For randomly-distributed O-vacancies, £ = 1/(5; however 
it is well known that O-vacancies cluster after annealing 
and that the chain fragments are typically much longer. 
In YBCO6.5, for example, the chains alternate between 
being completely filled and completely empty. In prac- 
tice, in any sample there will be a distribution of £, with 
some fraction n c h a in of these satisfying £ > 2^ c h a in- Since 
the magnitude of the low-field core contraction depends 
on the magnitude of the current circulating in the chain 
layer, we expect the low-field vortex core size to depend 
on n c hain- On the other hand, the crossover field B* de- 
pends primarily on the magnitude of the plane-chain cou- 
pling and should not be directly dependent on n c h a in^ 
Although data sets for 5 > 0.05 in Ref. [12] sample a 
limited set of magnetic field strengths, they appear con- 
sistent with this picture. 



V. CONCLUSIONS 

We have studied the vortex core structure within a 
simple proximity model for YBa2Cu307_^. We find that 
the current distribution around the vortex core is differ- 
ent in the chain and plane layers, and that the core is 
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elongated along the chain direction. There is a crossover 
in the magnetic field dependence of the core size at a field 
B*. The core size varies rapidly with magnetic field Bq 
for Bo < B* , and we have shown that B* is related to 
the energy scale of the small superconducting gap in the 
chain layer. Our calculations provide a natural explana- 
tion for the vortex core contraction measured in various 
/iSR experiments, and support the validity of the prox- 
imity model for YBa2Cu307_5. 
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